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i 
The outstanding features of this program are its flexibility in covering a 
very wide variety of problems and the multitude of program options 
available. The program is a combination of techniques and capabilities 
necessary in performing a complete aerodynamic analysis of supersonic and 
hypersonic shapes. These include: vehicle geometry preparation; computer 
graphics to check out the geometry; analysis techniques for defining vehicle 
component flow field effects; surface streamline computations; the shield- 
ing of one part of a vehicle by anorher; calculation of surface pressures 
using a great variety of pressure calculation methods including embedded 
H flow field effects; and computation of skin friction forces and wall 
temperature. 
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Although the program primarily uses local~slope pressure calculation methods 
that are mos. accurate at hypersonic speeds, its capabilities have been 
extended down into the supersonic speed range by the use of embedded flow 
field concepts. This permits the first order effects of component 
interference to be accounted for. 
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The program is written in FORTRAN for use on CDC or LBM types of computers. 


The program is documented in three volumes. Volume I is primarily a User's 
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FOREWORD 


Vhis report describes a compuler program developed at the Douplas 
Aircraft Division of the McDonnell Douglas Corporation, long Beach, 
California, The development of the Douglas Arbitrary-Body Aero- 
dynamic Computer Program was started in 1964 and greatly expanded 
in subsequent years under sponsorship of the Douglas Independent 
Research and Development Program (IRAD). From August 1966 to 
May 1967 the program development was continued under Air Force 
Contract No. F3361567-C-1008. ‘The product of thie work was the 
Mark il version of the program as released for use by government 
agencies in May 1967. Between 1967 and 1968 further Dougias IRAD 
work and another Air Force Contract (F33615-67-C1602) produced the 
Mark UJI Hypersonic Arbitrary-Body Program. The latest version of 
the program ay presented in this repurt is identified as the Mark IV 
Supersonic-Hypersonic Axbitrary-Body Computer Program and was 
prepared in the period of 1972-73 under Air Force Contract F33615- 
72-C-1675. This centract was administered by the Air Force Flight 
Dynamics Laboratory, Flight Mechanics Division, High Speed Aero 
Performance Branch. the Air Force Proiect Engineers for this study 
were Verle V. Bland Jr., and Captain Hugh Wilbanks, AFFDI./FXG. 


At the Douglas Aircraft Company, this work was conducted under the 
direction of Mr. Arvel E. Gentry as Principal Investigator. A number 
of major parts of the new program were prepared poy Mr. Douglas N. 
Smyth. Mr. Wayne R. Oliver's work in applying the various versions 
of this program to practical design problems contributed both in pro- 
gram design and in program validation. A number of other people 
contributed to the various phases of this work for which the authors 
are grateful. 


This report was submitted by the authors in November 1973. 


This technical report has besn reviewed and is approved. 


PHILIP P. 
Chief, Flight Mechanics Division 
Air Force Flight Dynamics Laboratory 
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ABSTRACT 


This report describes a digital computer prorram system that is 
capable of calculating the supersonic and hypersonic aerodynamic 
characteristics of complex arbitrary three-dimensiunal shapes. 
This program is identified as the Mark IV Supersonic-Hypersonic 
Arbitrary-Body Computer Program. This program is a com- 
plete reorganization and expansion of the old Mark III Hypersonic 
Arbitrary-Body Program. The Mark IV program has a number 
of new capabilities that extend its applicability down into the 
supersonic speed range. 


The outstanding features of this program are its flexibility in 
covering a very wide variety of problems and the multitude of 
program options available. The program is a combination of 
techniques and capabilities necessary in performing a com- 
plete aerodynamic analysis of supersonic and hypersonic shapes. 
These include: vehicle geometry preparation; computer graphics 
to check out the geometry; analysis techniques for defining 
vehicle component flow field effects; surface streamline computa - 
tions; the shielding of one part of a vehicle by another; calculation 
of surface pressures using a great variety of pressure calculation 
methods including embedded flow field effects; and computation 
of skin friction forces and wall temperature, 


Although the program primarily uses local-slope pressure calcu- 
lation methods that are most accurate at hypersonic speeds, its 
capabilities have been extended down into the supersonic speed 
range by the use of embedded flow field concepts. This permits 
the ficst order effects of component interference to be accounted 
for. 


The program is written in FORTRAN for use on CDC or IBM 
type compucers. 


The program is documented in three volumes. Volume I is pri- 
marily a User's Manual, Volume II gives the Prograrn Formulation, 
and Volume III contains the Program Listings. 
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SECTION I 
SUMMARY 


‘This volume contains the source language listings of the Mark IV 
Supersonic-Hypersonic Arbitrary-Body Program (Mod 0 Version). 
The program as shown in this listing will operate on CDC 6500, 
6600, and CYBER 74 computers. With a small converter program, 
the Mark IV program can be converted for operation on IBM 360 and 
370 types of computers. This converter program is included with 
the listings. 


The actual tapes released for this program may contain minor 
differences from these listings. Therefore, these listings should 

be used for reference and communication purposes only. Any 
problems, bugs, etc., encountered with this pyogram should be 
brought to the attention of the authors as soon as possible. Also, 
because of the complexity and size of the program, users should 
consult with the authors before making extensive program additions 
or modifications. This will sometimes save time since the proposed 
changes may have already been accomplished by someone else, 
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SECTION II 


DATA STORAGE TECHNIQUES 


Each major part of the Mark IV program makes u&e of data that has been 
generated and stored by some other program component. The manner in 
which these data are stured and retricved is important both to the user 


and to the programmer seeking to ada some new capability to the program. 


The principal concern for the user is that he know what is stored, and 
what sequence number, pointer, etc., might be required as subsequent 
input data to retrieve it. The primary storage units are tabulated below. 


Unit 
Number Name Usage 
I TAPEIN Copy of all input data. 
Pa TAPEB Scratch unit (used in SHIELD and GENCUT). 
3 TAPEA Negative elements from shielding analysis. 
4 - Geometry data storage (mass storage). 
5 TAPEIN System input data stream (copied to unit 1 
in routine MONITR, TAPEIN set = 1). 

6 TAPEOT = System output print unit. 

7 TAPEG System punch output unit. 
8 - Element data geometry data storage (Type 3 cards) 
9 - Force data storage (mass storage) 

10 - Flow field data storage (mais storage) 

1 TAPEC Scratch unit (used in AUXIL.I) 

12 TAPEE Geometry data for interference usage. 

13 TAPEF Geometry data plus local flow field data. 


14 TAPED Scratch. 


Units 4, 9, and 10 listed above are identified as being ''mass storage" 
units. That is, the units are used to store data along with the necessary 
tags or pointers so that a given item can be retrieved directly without 
using a sequential read. On the CDC machine this process is called 
"mass storage", On the IBM machines it is identified as "direct access" 
storage. Although the two machines use different type of FORTRAN 
statements to write and read the data, the general procedure is still the 
game. A pointer or tag must be set or calculated by the program for 
each FORTRAN write statement. This pointer or tag is the sequential 
humber of the record in ths data set. The use of the same tag number in 
a FORTRAN read will result in the proper retrieval of the record. The 
apecific data stored on the mass storage units 4, 9, and 10 are described 
below. 
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Unit 4. Geometry Data 


Unit 4 is the storage unit for geoinetry data generated in routine GEOM. 
The unit normally contains both the input element corner points plus all 
the characteristics of the plane quadrilateral. When it is requested ty 
do so, the program will also store on unit 4 the surface property data 
calculated in the FORCE routine. Unit 4 has a master table of contents 
followed by the data fo. each element. The unit is set up and initiaiized 
in routine GEOM. 


Record Word 


Nurnber Number Description 
1 1 Configuration title (CONFIG from the Geometry 


Contrvul Card) 


NEXT, pointer to next empty recore. 

NP, number of PANELS currently stored. 

NPMAX, maximum number of panels provided 
for in the table of contents. 

4 NREM, number of record spaces provided for use 

on each element (set = 3 in GEOM). 


wn 


3 1 Pointer to next record after geometry data storage 
is complete. 


Not used. 


5 Panel #1 table array. 

Panel number (ISTAT3 = 1). 

Panel name (PANEL from the Panel Identification 
Card) 

Record number of first element (ISTART). 
Number of elements in the panel (L). 

Panel orientation (IORN) 

Symmetry flag (SYMFCT} 

Number of columns of elements (N) 

Number of rows in each column (I) 

Panel number adjacent to side | of panel (NADJ1) 
10 Panel number adjacent to side 2 of panel (NADJZ2) 
11 Panel number adjacent to side 3 of panel (NADJ3) 
12 Panel number adjacent to side 4 of panel (NADJ4) 


Blank 
Blank 
Blank 
Blank 


Panel #2 table array 

(same format as record number 5) 
11 Blank 

lé Blank 

13 Blank 

14 Blank 


15 Fanel #3 table array 
(same format as record number 5) 
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Record Word 





Number Number Description 

(NP*5) Panel #NP (last panel lvaded) 

(NP*5)+1 Blank 

(NP*5 )+2 Blank 

(NP*«5)+3 Blank 

(NP«5)+4 Blank 

INPMAX*5)+1 Last panel for which space is provided 

(NP MAX*5)+2 (NPMAX may be = NP if desired) 

(NPMAX*5)+43 

(NPMAX*5)+4 

{((NPMAX#1)*5] Element #1, First available record tor element 
data. The array ELEM is written here in GEOM, 

{(NPMAX+1)*5+1] Pointer array (25 words) to where surface output 

data is written by FORCE for input side of vehicle. 
First word is for lst a-8, 2nd for 2nd a-8, etc. 

[((NPMAX#1)*5+2] Same as above for reflected side of vehicle. 


The data written on Unit 4 for each element are contained in the element 
data array, ELEM, in subroutine GEOM. As indicated above, the 25 
word array EJEM for the first element is written on Unit 4 at record 
number [(NPMAX+1)*5]. In routine GEOM the next two records, 
[(NPMAX4+])*5+1] and [(NPMAX+1)*5+2], are filled with 0.0 values. These 
records will be filled with live pointer data in subroutine FORCE as indi- 
cated above. These pointers will give the location of the calculated force 
data. We thus have three records prr element on Unit 4. 


The element data array, EILEM, contains the following information: 


ELEM(1) 
ELEM(2) 
ELEM(3) 
ELEM(4) 
ELEM(5) 
ELEM(6) 
ELEM(7) 
ELEM(8) 
ELEM(9) 
ELEM(10) 
ELEM(11) 
ELEM(12) 
ELEM(13) 
ELEM(1\4) 
ELEM(15) 
ELEM(16) 
ELEM(17) 
E LEM(18) 


L. Element number 

N Column number 

I Row number 

NX X-component of unit normal 

NY Y-component of unit normal 

NZ Z-component of unit normal 
XCENT X-coordinate of clement centroid 
YCENT Y-coordinate of element centroid 
ZCENT 2Z-coordinate of element centroid 
AREA Element area 
XIN(1) X-coordinates of input elernent 
XIN(2) corner points 
XIN(3) 
XIN(4) 
YIN(1) Y-coordinates of input element 
YIN(2) corner points 
YIN(3) 
YIN(4) 
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ELEM(19) = ZIN(1) 2Z-coordinates of input element 
ELEM(20) = ZIN(2) corner points 
ELEM(2]) = ZIN(3) 
ELEM(22) = ZIN(4) a 
ELEM(23) = 0.0 
ELEM(24) = 0.0 
7 0.0 


f 
ELEM(25) 
tO 


When directed by the input parameter ISAVE subroutine FORCE will 
write certain results from the pressure and force calculations out on 
Unit 4. These surface property data are in the form of single records 
containing the array El in FORCE. Since up to 20 @-8 conditions can 
be analyzed on a pass into FORCE, pointers are required on Unit 4 to 
indicate where each El record is written, These pointers are written 
in the vacant record locations provided righi behind the original element i 
data record, ELEM, for each element. The surface property data 1 
contained in the array El in FORCE are as follows. 





El(l) = IG4 Pointer to original element data on Unit 4 
El(2) - E(1) Element number : 
E1(3) = 0.0 
El({4)) = 0.0 
E1(5) = cP Element pressure coefficient 
E1(6) = DELTA Element impact angle, . ¥ 
El(7) = Sx X-component of surface unit velocity vector | 
E1(8) = SY Y-component of surface unit velocity vector a 
El1(9) = SZ Z-component of surface unit velocity vector = 
E1(10) = BS(2)/PFS Local pressure divided by freestream Bro! 
pressure ote : a 
El(il) = BS(3)/TFS Local temperature divided by freestream ES 
temperature | @ 
El(1?) = DELCA AC, | a 
E1(13) = DELCY ACy q 
E1(14) = DELCN AC, | We 
E1l(i5) = DELCLL ACg pas 
E1(16) = DELCLM AC, a 
: El(17) = DELCLN AC, an 
| E1(18) = 0.0 § 
| E1(19) = 0.0 | 
E1(20) = 0.0 ‘ : 
E1l(21) = 0.0 i 
El(22) = 0.0 ; 3 
E1(23) = 0.0 A 
‘ E1(24) = 0.0 : 
E.1(25) = 0.0 
The total number of records that can be stored on Unit 4 is controlled 4 
by the pointer array size provided for in the MAIN program. In the { 


initial release c. the program this was set at 3000 records. A check is 
inade in routine FORCE to prevent the storage of more than this number of 
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records. Thus, the total number of elements that can be analyzed at 
one time depends upon whether or not surface property data are saved. 


The total number or record spaces required on Unit 4 is given by the 
following relationship. 


L = 4+ NPMAX*5 + L *3 4+ L 
max 


total save 


where 
NPMAX 


Maximum number of panels provided 
for on the unit (input on the Geometry 
Control card). 


L = The total number of elements for all 
imax 
panels. 
L = The tota! number of elements for 
save 


which local surface property data 
will be saved in routine FORCE. 
Must include the total element count 
for both sides of the vehicle if it is 
yawed and the symmetry reflection 
process is used. 


The size of the pointer storage array for Unit 4 in the MAIN program 
should be increased from 3000 up to 6000 to 10000 if additional core 
space is available on the user's computer installation. 


Unit 9. Force Data Storage 


Unit 9 is used to store component force data generated in the inviscid 
routine PRES, and in the viscous routine VISCUS. Unit 9 contains a 
table of contents followed by the data for each vehicle component. 


Record Word 
Number Number Description 


l l NTOT, total number of components. 
2 ISUM, total number of summations. 
3 NEXT, next empty record space. 

4 TITLE9(1) = LOHFORCE DATA for CDC 
= 4HFORC for IBM 
5 TITLE93(2) = 1OHSAVE UNIT for CDC 
= 4HE DA for IBM 
6 TITLE9(3) = blank for CDC 
= 4HTA S for IBM 
7 TITLE9(4) = blank for CDC 
= 4HAVE for IBM 
8 TITLE9(5) = blank for CDC 


4HUNIT for IBM 
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Record Word 
Number Number Description 


2 1 Identification flag for component data (=0.0 for - 
force data only, = 1.0 if derivatives are included. 
2-41 LCOM, pointers to location of component data 
(40 total). 


3 1 Identification flag for summation data (= 0.0 for 
force data only, = 1.0 if derivatives are included). 
2-41 LSUM, pointers to location of sumrnation data 


(40 total). 








4 1 Mach, freestream Mach number. 
2 V, freestream velocity. 
3 RENO, freestream Reynolds number per foot. 
4 ALT, altitude. 
5 PSTAG, stagnation pressure in atmospheres. 
6 TSTAG, stagnation temperature. 
7 GTYPE(1), gas type name. 
8 GTYPE(2), gas type name. 
9 SREF, reference wing area. 
10 SPAN, reference span length. 
11 MAC, reference mean aerodynamic chord. 
12 XCG, X-location of cg. 
13 YCG, Y-location of cg. 
14 ZCG, Z-location of cg. 
5-9 Not used. 
10 Force data table of contents for lst component 
(record number is = LCOM(1) = 10). 
1 Component number = 1. 
2-11 IPANL array, panel numbers used to make up 
this component. 
11 1 NAB, number of a@-f sets. | 
2 Value of @ for lst @-8 set. q 
3 Value of a for 2nd a-B set. j 
etc. 
21 Value of B for lst @-£ set. { 
22 Value of B for 2nd a@-B set. ~ 
etc. 
41 2 i 
12 Force data for lst component at lst a-£. ; 
The array F in routine PRES. : 
1 F(l) = @ ; 
3 F(3) = Cy, i 
4 F(4) = Ca { 
5 Fi5) = Cy 
6 F(6) = 
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Record Wore 
Number Number Description 


7 F(7) = B 
8 F(8) = Cy/Cp 
9 F(9) = m 
10 F(10)= Cy 
ll F(Ll) = Cy 
13 1-11 Force derivatives for 1st component at ‘st a- £6. 
Not used. 
14 1-11 Force data for lst component at 2nd a- 8. 
Array F. 
15 1-11 Force derivatives for lst component at 2nd a- B. 
Not used, 
16 1-11 Force data for lst component at 3rd a-£f. 
Array F. 
17 1-11 Force derivatives for lst component at 3rd a- @. 
Not used. 


etc. 


The above storage process starting with record 10 is repeated for each 
vehicle component. The starting record number for each component is 
placed in the LCOM array that is stored in record number 2. The 
pointers to summation data sets are placed in LSUM and stored in 
record number 3. Normally, record number 4 containing the Mach 
number and other data is based on the last component placed on the unit. 


Unit 10. Flow Fiel¢ Data Storage 


Unit 10 is the most complex of the mass storage (direct access) units 
used by the prcgram. Although Unit 10 is called the Flow Field unit it 
is used to store more than just flow field data. It is also used to store 
surface property data, streamline data, and skin friction computation 
results. The unit contains a number of directories and tables that hold 
not only basic identification information, but pointers that provide 
storage location record numbers pointing to lower level directories 
and the various types of data. 


The general format of Unit 10 is shown in Figure 1. Starting with the 
Master Director at the top we see that it contains a set of 5 pointers in 
the array IMTAS8. These point to what are called major data sets. For 
example, each set could be for a completely different vehicle or for a 
different Mach number for the same vehicle. Each of these sets has a 
Set Directory containing the Mach number and a pointer array, LOAB, 
pointing to a different lower level director for each a@-B combination. 

In turn, each a- director has pointers, LORG, that point to the Region 
Directories. The Region Directories contain some basic flag ana inter- 
polation parameters plus a set of pointers used to reach and retrieve 
data in the still lower level sub-regions and secondary flow field tables. 
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Record 
Number Description 
Master Directory 
1 TITLEM = Unit title (10 words) 


2 IMTAB(1) = Number of separate data sets stored on 
the unit (NDSET) 


1 
i 
! The basic format for each of the directory tables is listed below. 
{ 
| 
j 


IMTAB(2) = Pointer to next empty record (LCNEXT) 

: IMTAB(3) = Empty 
IMTAB(4) = Empty 
IMTAB(5) = Pointer to start of lst set = LOSET(1) 
IMTAB(6) = Pointer to start of 2nd set = LOSET(2) 
IMTAB(7) = Pointer to start of 3rd set = LOSET(3) . 
IMTAB(8) = Pointer to start of 4th set = LOSET(4) | 
IMTAB(9) = Pointer to start of 5th set = LOSET(5) 7 

3 Empty 


| Set Directory (First set is at record number 4) 
| LOSET TITLES, MACH, NAB (number of @-f sets) 

| LOSET+1 WLOAB array pointers to a-@ sets 

| LOSET+2 Empty 

i 


a-§ Directory 


LOAB TITLEA, ALPHA, BETA, NREG (number of regions) 
LOAB+t+l1 LORG array pointers to region directories 
LOAB+2 Empty 


Region Directory 


1 uniform flow field data 
2 non-uniform data 
ITFLAG Not used 


LORG+tl XO, YO, ZO, PSIO, etc. q 


LORG+2 PHI or YPA (not active) 
LORG#+3 PHI or YPA (not active) 


LORGt4 Empty 


LORG+5 Pointers for sub-region #1 i. 
LOFF, IFC, LOCD, IFD, LOSF PB 
LORG+6 =Pointers for sub-region #2 ,. 
LOFF, IFC, LOCD, IFD, LOSF Py 

etc. etc. 
LORG+(4+NSRG) Pointers for last sub-region (NSREG) 


LORG TITLER 
NSREG = Number of sub-regions 
IDTYP(1) = Data type flag 
= 1 flow field data 
= 2 surface property data 
= 3 streamline data 
= 4 shock shape data (not used at present time) 
IDTYP(2) = Uniform - non-uniform flag 
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SECTION Iil 


IBM GO-DECK SETUP 


A suggested set of DD cards for operation of the Mark IV program 
on IBM 360 and 370 computers is given below. 


// EXEC FORTGO, REGION=2 40K 


//STEPLIB DD  ODSNAME=SYS3.C2BC, DISP=(SHR, PASS), UNIT=3330, 


// VOLUME=SER=CSLB07 

//GOs-FTO1F001 OLD DSN=8ST01, UNI T=SYSDA, SPACE=C(TRK,C5,1)), 

// DCH=CRECFM=F, LRECL=80, BLKSIZE=80) 

//GO.FTO2F001 DD DSN=§£T02, UNI T=SYSDA, SPACE=C(TRK,C10,5)), 
// DC B= CBUFNO=1, BLKSIZE=6447, RECFM=VBS) 

//GO.FTO3F001 OD DSN=§FT03, UNIT=SYSDA, SPACE=C(TRK,(€10,5)), 
fi DC B=(BUFNO=1, BLKSIZE=6447, RECFM=VBS) 

//GO.FTO4F001 DD OSN=§FTO4, UNI T=SYSDA, SPACE=(100, €3000)), 
// DC B= CBUFNO=1 , DSORG= DA) 

//GO.FTN7F001 OD SYSOUT=B 

//GO.FTOBFO01 DC DSN=§FTO8, UNIT=SYSDA, SPACE=C(TRK,(C25,5)), 
// DCB=CRECFM=F, LRECL=80, BLKSIZE= 80) 

//GO.FTUSFO001 DD DSN=EFT09, UNIT=SYSDA,SPACE=(164,(500)), 

// DC B= ( GUFINO=1, DSORG=DA) 

//GO.FT10F001 DD OSN=&FT19, UNI T=SYSDA, SPACE=(100, €2000)), 
// DCB=( BUFNN=1, DSURG= DA) 

//GO.FT11F901 DD DSi EFT11, UNI T=SYSDA, SPACE=CTRK,(10,5)), 


// DCB=(BUFNO=1, BLKSIZE=6447, RECFM=VBS) 
//GO.FT12F001 DD DSH=EFT12, UNIT=SYSDA,SPACE=CTRK,C10,5)), 
// DCB=( BUFNO=1, BLKS1ZE=6447, RECFM=VBS) 
//GO.FT13FQ01 DD DSN=&FT13, UNI T=SYSDA, SPACE=C(TRK,(C10,5)), 
// DCB=CBUFNO=1, BLKSIZE=6447, RECFM=VBS) 
/7/GO.FT14F001 DD DSN=§FT14, UNIT=SYSDA, SPACE=C(TRK,(€10,5)), 
/! DCB=( BUFNO=1, BLKSIZE=6447, RECFM=VBS) 


//GO.SYSUOUMP DD SYSOUT=A,SPACL=CTRK,(100),RLSE) 
//GO.SYSIN DD * 


ll 
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SECTION IV 


OVERLAY STRUCTURE 


Because of the large size of this program it is necessary to use the over- 
lay features of FORTRAN for normal operation. The overlay structure 
for the CDC version of the program (CDC SCOPE 3.3 overlay) is shown 
below. With thia structure tho program requires approximately 114,100 
octal words on the CDC machine (CYBER 74, FTN V3.0-P308 Opt=1). 


CDC Overlay 


4 OVERLAY (MARK4 , 0,0) 
MAIN 
HEADER 
MONITR 
ARSIN 
ARCOS 

OVERLAY (MARK4, 1,0) 
GEOM 
PUNCH 
OVERLAY (MARK4, 1,1) 
IELE 
OVERLAY (MARK4 , 1,2) 
ELLIP 
OVERLAY (MARK4 , 1, 3) 
CUBIC 
OVERLAY (MARK4 ,1,4) 
AIRCFT 
SURF 
NACEL 
FUSE 
STATS3 
OVERLAY (MARK4, 2 ,0) 
AERO 
COMPR 
EXPAND 
CONE 
NEWTPM 
ATMOS 
SOLVIT 
GETT 
SAVE 
ROWFM1 
OVERLAY (MARK4 , 2,1) 
FLOW 
MERID 
NCONE 
CONEA 
MOUT 
CADA 
FFSPEC 
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FFAIR 
FFSURF 
FFINPT 
FFBODY 
SOSE 
OVERLAY (MARK4 , 2,2) 
SHIELD 
OVERLAY (MARK4 , 2, 3) 
PRES 
FORCE 
SHKEXP 
BLUNT 
CPINPT 
FLNTRP 
VALU2 


ACONE 
OVERLAY (MARK4 , 2, 4) 
VISCUS 
FORSF 
SKINFR 
TEMP 


ROMU 3 
OLY 
piock DATA 

| 


a 


CFINPT 


SPLINE 
OVERLAY (MARK4 ,2,5) 
SPEC 
SUM 
OVERLAY (MARK4 ,2,6) 
STREAM 
SENTRP 
VALUE 
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OVERLAY (MARK4 , 3,0) 
GRAPH 

OVERLAY (MARK , 4 ,0) 
AUXILT 
GENCUT 
OUTD 


IBM Overlay 


The overlay for the IBM version of the program is shown below. 
With this structure the program requires approximately 240K bytes 
(FORTRAN IV H Opt-2, Level 20.1, OS/MVT Release 201). 


INSERT MAIN,HEADER, EROR,MONITR 
INSERT EXEC, TAP, TAGS 
OVERLAY ALPHA 
INSERT GEOM 
INSERT GFLAG 
OVERLAY BETA 
INSERT IELE,ELLIP,CUBILC,PUNCH 
OVERLAY BETA 
INSERT ALRCFT, SURF, FUSE ,NACEL STATS 3 
INSERT ACFT,ARFOIL, FUSD,NACELD 
OVERLAY ALPHA 
INSERT AERO,ATMOS 
INSERT FLIGHT,GASD, ABDATA, REF, INTERF , FSBS 
INSERT SHKEXP ,COMPR, EXPAND ,NEWTPM, CONE , BLUNT ,GDATA, FARRAY ,MDATA 
INSERT SOLVIT ,ROWFM1 , SAVE, GETT 
OVERLAY BETA 
INSERT FLOW, FFINPT , FFBODY ,FFAIRF ,FFSPEC, FFSURF 
INSERT SOSE,NCONE,CONEA,MERID,MOUT 
INSERT FFIELD,SOSEIN,SOSEOT,CADA 
OVERLAY BETA 
INSERT SHIELD 
OVERLAY BETA 
INSERT PRES, FORCE, ACONE 
INSERT FLNTRP,CPINPT , SFNTR2, SURFN2, VALU2 
OVERLAY BETA 
INSERT VISCUS , FORSF, TEMP ,QC, POLY ,PROP, ROMU, TEMPQC ,SKIND,C1 
OVERLAY BAMA 
INSERT SKINFR 
OVERLAY GAMA 
INSERT INTEG, RUNKUT ,GRADNT , LGRNGE, FUNCT, INTL, INT2,SMTHNA, SPLINE 
INSERT ROOT, CURVFT, SIMPS}. 
OVERLAY DELTA 
INSERT PRECAL 
OVERLAY DELTA 
INSERT LAMNAR 
OVERLAY DELTA 
INSERT TURBLN,PROFIL 
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OVERLAY GAMA 
INSERT CFINPT,SFNTR3,VALU3 
OVERLAY BETA 
INSERT SPEC, SUM 
OVERLAY BETA 
INSERT STREAM 
INSERT SENTRP, VALUE, SURFNT 
OVERLAY ALPHA 
INSERT GRAPH 
OVERLAY ALPHA 
INSERT AUXILI 
INSERT GENCUT,OUTD 


1 
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The listing as presented hero is for the CDC version of the program. This 
program has been run on CDC 6500, 6600, and CYBER 74 models. The 
program is written in FORTRAN for the CDC FTN compiler (opt=]) and 
hes been run under the SCOPE 3.3 and 3,4 operating svstems. 


In this listing all cards that are peculiar to the CDC versionof FORTRAN 
are identified by a C in card column 80. All carda that are peculiar to 
the 1BM FORTRAN IV compiler are identified by an I in card column 80 
anda C in card column 1. Jn other words, the code for both CDC and 
BM machines is in the deck but the IBM cards are made inactive by con- 
, verting them to comment statements (C in card column 1). Since all of 
the machine dependent cards are identified by an I or C in card column 
80 it is a simple matter to convert the deck from one version to the other 
| with a small conversion program. When converting from CDC to IBM 
code this conversion program reads and copies. ch card toa storage 
; unit. If a card has a C in card column 80, then a C is written into card 
j column 1 to make the CDC peculiar card inactive. If a card has an [in 
| card column 80, then the C is removed from card column 1 and the card 
image wriiten to the storage unit as an active FORTRAN statement. The 
conversion from IBM back ta CDC is made in a similar manner. 


The conversion program to convert the deck from CDC to IBM FORTRAN 
is listed below (for use on an IBM machine). 


DIMENSION DATA(22) 
DATA CB,CC,CI/1H ,1HC,1HI1/ 
REWIND 19 
DO 100 I=1,20000 
KEAD (5,20,END=300) DATA 

20 FORMAT (J.A1,19A4,1A2,1A1) 
IF (DATA(2Z) .EQ. CI) DATA(1) = CB 
IF (DATA(22) .EQ. CC) DATA(1) = CC 
WRITE (19,20) DATA 

100 CONTINUE 

300 STOP 
END 
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! 
APPENDIX A 
PROGRAM LISTINGS 
This part of the report contains the source card listings for the Mark IV 
Mod 0 Version of the Supersonic-Hypersonic Arbitrary-Body Program. 


This program places the new deck with IBM cards made active, and CDC 
cards inactive, onto unit 19, The references to unit 19 above can be 
changed to unit 7 to punch the deck out. 
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Routine 


MAIN 
MONITR 
HEADER 
EROR 
ARSIN 
ARCOS 
GEOM 
PUNCH 
TELE 
ELLIP 
CUBIC 
AIRCFT 
SURF 
NACEL 
FUSE 
STATS3 
AERO 
COMPR 
EXPAND 
CONE 
NEWTPM 
ATMOS 
SOLVIT 
GETT 
SAVE 
ROWFMI1 
FLOW 
ME RID 
NCONE 
CONEA 
MOUT 
CADA 
FFSPEC 
FFAIRF 
FFSURF 
FFINPT 
FFBODY 
SOSE 
SHIELD 
PRES 
FORCE 
SHKEXP 
BLUNT 
CPINPT 
FLNTRP 
VALUE 
SFNTR2 


INDEX TO SUBROUTINE LISTINGS 


Name 


Main Program (Executive Program) 
Input Data Read and Frint Routine 

Page header subroutine 

IBM ABEND dump routine 

Arcsine conversion routine (CDC) 
Arccosine conversion routine (CDC) 
Geometry Data Preparation and Control 
Element data write subrou:ine 

Element data transfer routine 
Elliptical cross-section routine 
Parametric cubic routine 

Aircraft goometry executive program 
General airfoil surface routine 

Nacelle geometry generation routine 
Fuselage geometry generation routine 
Dummy zero area element routine 
Aerodynamic analysis executive routine 
Compression routine 

Prandtl-Msyer expansion routine 
Conical flow routine 

Blunt body Newtonian + Prandtl- Meyer 
Atmosphere subprogram 

Matrix operation routine 

Data retrieval routine for SOLVIT 

Data save routine for SOLVIT 

Matrix formation for surface spline 
Executive routine for flow field analysis 
Meridian cutting plane routine 
Improved cone flow routine 

Cone at angle of attack routine 
Meridian plane output data organization 
Circular arc data analysis routine 
Simple flow field data input routine 
Dummy routine for future use 

Surface data transfer routine 

Flow field data hand load routine 

Body flow field data generation routine 
Second-order shock-expansion routine 
Shielding routine 


Executive routine for pressure calculations 


Inviscid pressure and force calculation 
Element strip shock expansion routine 
Blunt-body viscous routine 

CP surface interpolation routine 

Flow field interpolation routine 

Spline interpolation routine 

Surface data interpolation routine 
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Deck 


MAIN 
MONI 
MEAD 
EROR 
ARSI 
ARCO 
GEOM 
PUNC 
IELE 
ELLI 
CUBI 
AIRC 
SURF 
NACE 
FUSE 
STAT 
AERO 
SOMP 
EXPA 
CONE 
NEWT 
A‘CTMO 
SOLV 
GETT 
SAVE 
ROWF 
FLOW 
MERI 
NCON 
CONE 
MOUT 
CADA 
FFSP 
FFAI 
FFSU 


“FFIN 


FFBO 
SOSE 
SHIE 
PRES 
FORG 
SHKE 
BLUN 
C PIN 
FLNT 
VAL2 
SFNe 


Page 


12 
13 
14 
15 
16 
29 
31 
33 
41 
49 
62 
73 
77 
82 
83 
88 


97 
102 
108 
112 
125 
126 
127 
130 
140 
155 
158 
163 
172 
175 
179 
180 
188 
196 
206 
223 
250 
261 
284 
289 
293 
304 
326 
332 








Routine 


ERF 
ERRF 
ACONE 
VISCUS 
FORSF 
SKINFR 
TEMP 
ROMU 
Qc 
POLY 
BLOCK 
CFINPT 
INTEG 
PRECAL 
SFNTR3 
VALU3 
LAMNAR 
PROFIL 
TURBLN 
RUNKUT 
FUNCT 
INT1 
INT2 
LGRNGE 
ROOT 
CURVFT 
GRADNT 
SMTHNA 
SIMPS1 
SPLINE 
SPEC 
SUM 
STREAM 
SFNTRP 
VALUE 
GRAPH 
AUXILI 
GENCUT 
OUTD 
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Name 


Error function routine (CDC) 

Error function calculation (CDC) 

Cone at angle of attack 

Viscous executive routine 

Skin friction force routine 

Mark IIL skin friction calculation routine 
Surface temperature calculation routine 
Density-viscosity routine 

Heating at given wall temperature routine 
N-th order polynomial routine 

Block data routine for skin friction 
Sxin friction data interpolation control 
Executive routine for integral boundary lay. 
Preliminary calculation routine 

Skin friction interpolation routine 

Skin friction interpolation 

Integral method laminar boundary layer 
Velocity profile calculation 

Turbulent boundary laver routine 
Integration of turbulent equations 
Integral boundary layer functions 
Integral method function 1 

Integral method function 2 

Lagrangian interpolation 

Root finding routine ‘ 
Curvesefit routine 

Curve gradient routine 

Curve smoothing routine 

Simpsons rule integration routine 
Spline curve fit routine 

Special function routine 

Force data summation routine 
Streamline calculation routine 

Surface data interpolation 

Spline interpolation routine 

Graphics executive routine 

Auxiliary methods executive routine 
General cutting plane rovtine 

Cutting plane output data -:ollection 


Deck 


ERF 
ERRF 
ACON 
visc 
FORS 
SKIN 
TEMP 
ROMU 
QC 
POLY 
BLKD 
CFIN 
INTE 
PREC 
SFN3 
VAL3 
LAMN 
PROF 
TURB 
RUNK 
FUNC 
INT 1 
INT2 
LGRN 
ROCT 
CURV 
GRAD 
SMTH 
SIMP 
SPLI 
SPEC 
SUM 
STRE 
SFNT 
VALU 
CRAP 
AUXI 
GENC 
OUTD 


Page 


342 
343 
344 
350 
360 
370 
381 
389 
394 
398 
399 
402 
414 
429 
435 
445 
452 
464 
469 
472 
477 
479 
480 
482 
484 
485 
486 
487 
488 
490 
491 
492 
500 
516 
526 
$32 
533 
535 
551 
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